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semiminor axis of ellipse 

body bluntness parameter for conic section ((b/a) 2 for ellipse) defined 
by y 2 = 2R b x - E^x 2 (ref. l) 

shock bluntness parameter for conic section defined by y 2 = 2R x-B x 2 
(ref. l) s s 

semimajor axis of ellipse or span measured from shock apex 

constant in stream function (eq. (Cl)) 

chord or body length 

length along chord normalized by c 

scale factor (eqs. (D7), (D8)) 

stagnation point mass-flow function (eq. (3)) 

constants in velocity distribution functions (eqs. (39)- (^2)) 

Mach numb er 

location of control surface (sketches (e) and (f)) 

correlation exponent (eq. (l5)) 

pressure on body 

streamline radius behind shock 

wing section or body radius at stagnation point 
shock radius at apex 

edge or corner radius — - 

base radius of conical body (fig. 15) 

body radius of section between tangent points (sketch (j)) 
distance on body surface from stagnation point 
distance from stagnation point normalized by s* 
velocity component in 5»n coordinate system (sketch (h)) 
velocity 


iii 



V 


local stream velocity normalized by V 

Y derivative of V with respect to s 

AV correction term to velocity distribution function (table 2) 

v velocity component in £,r> coordinate system (sketch (h)) 

v , dimensionless Newtonian velocity 

Av' ,Av" nonlinear components of v (sketch (d)) 
v base (linear) component of v (sketch (d)) 

-D 

v derivative of v with respect to s 

X,Y coordinate system defined in figure 10 

X g ,Y s shock apex coordinates, figure 10 

x,y coordinate system with origin at shock apex . 

Y reference ordinate, figure 10 
m 

y vertical distance at shock between shock apex and stagnation 

streamline (sketches (b) and (i)) 

y vertical distance at body between shock apex and stagnation 'point 

(sketches (b) and (i)) .. .. . . - 

y* vertical distance between shock apex and sonic point on body 

(sketch (a)) 

a angle of attack 

Y specific heat ratio 

A shock standoff distance 

6 stream deflection angle behind shock 

e angle of tangency point of edge radius with respect to chord plane 

normal (sketch (c) ) . .. . ...... - 

r, vorticity . .. ... 

0 shock surface . inclination with respect to plane normal to.free- 

stream direction; in degrees , except where defined otherwise 

£,n coordinate system with origin of stagnation point (sketch (h)) 

P 

iv 


gas density 


density of gas behind shock inclined at angle 0 

average mass flow per unit area passing a specified control surface 

body surface inclination with respect to plane normal to free-stream 
direction; in degrees, except where noted 

stream function (eq. (Cl)) 

stagnation stream impingement angle (sketch (h)) 
correction term to trial shock apex location 

Subscripts 

cylinder 

lower 

normal to control surface 
plate 

at shock location where streamline has the radius of 
at stagnation point or along stagnation streamline 
upper 

at location on body in free stream in direction from shock apex 

free stream 

behind normal shock 


Superscripts 

at sonic point location 

tangent point location on body (sketch ( j ) ) 


A METHOD FOR PREDICTING SHOCK SHAPES AND PRESSURE DISTRIBUTIONS 


ON TWO-DIMENSIONAL AIRFOILS AT LARGE ANGLES OF ATTACK 
George E. Kaattari 
Ames Research Center 


SUMMARY 


A method is presented for determining shock envelopes and pressure distri- 
butions for two-dimensional airfoils at angles of attack sufficiently large to 
cause shock detachment and subsonic flow over the windward surface of the air- 
foil. Correlation functions obtained from exact solutions are used to relate 
the shock standoff distance at the stagnation and sonic points of the body 
through a suitable choice for the shock shape. The necessary correlation func- 
tions were obtained from perfect gas solutions but may be extended to any gas 
flow for which the normal shock-density ratio can be specified. 


INTRODUCTION 


There is current interest in the aerodynamics of planetary atmospheric 
entry vehicles. Such vehicles are blunted to ensure a detached shock envelope 
in order to minimize vehicle heating rates during the flight trajectory. 

Most theoretical and experimental investigations have been devoted to the 
prediction and measurement of shock envelopes and pressure distributions about 
axisymmetric shapes. Some attention has also been given to a limited class of 
two-dimensional cylindrical bodies (ref. l) and to two-dimensional asymmetric 
flow about flat plates at large angles of attack (ref. 2). However, solutions 
for asymmetric two-dimensional shapes are not readily available. 

This investigation was undertaken to develop an approximate method for 
estimating the shock envelope and pressure distribution over the windward sur- 
face of a large variety of two-dimensional profiles at angles of attack suf- 
ficiently large to result in shock detachment. The method employs a technique 
similar to that described in reference 3 but utilizes correlations based on 
available exact solutions for two-dimensional flows. 


ANALYSIS 


The analysis will be developed in four parts. First, the effect of body 
shape on the velocity gradients and shock standoff distance along the stagna- 
tion streamline is considered. Second, basic relationships between shock and 
body geometry in symmetric flow are derived and applied to symmetric flow shock 
solutions. Third, the symmetric flow relationships are modified and extended 



to asymmetric flow shock solutions and demonstrated with a numerical example. 
Finally, the method for estimating the pressure distribution is developed and 
demonstrated by a numerical example. 

Computing forms and necessary figures are provided to organize and expe- 
dite both the shock and pressure distribution calculations. 


Stagnation Point Velocity Gradients 


The effect of body radius on the velocity gradients and the shock standoff 
distance along the stagnation streamline is established as follows: Traugott 

(ref. 4) showed that the dimensionless velocity gradient at the stagnation 
point on a blunt body is a function only of the free-stream flow properties , 

Mj and y: 


A (dV/ds) 

~ (1 + VV ’ f(Ml,Y) 


(1) 


Assuming that this correlation is applicable to all shock solutions, any par- 
ticular solution can be used to determine f(Mi,y). Accordingly, the solutions 
in reference 1 for shock-standoff distances, A ^/R^, and pressure distributions 
on circular cylinders were used to extract numerical values for the velocity 
gradient in the vicinity of the stagnation point for various specific-heat 
ratios, y, and free-stream Mach numbers. Substituting this gradient into equa- 
tion (l), which can be rearranged into a convenient form, yields numerical 


values for f(Mi,y): 


f(Mi,y) = 


jst_ \ 
Rb Vi 


dV 
. ds 


st 



(2) 


For purposes of the flow continuity analysis described in appendix A, 
equation (2) is modified to obtain the stagnation point mass-flow gradient: 
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<a(p st v/PiV!) 

P st 
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d(s/A ) 
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Pi 
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f(M ls y) 


= G 


(3) 


The function G is plotted as a function of the normal shock-density ratio, 
P 2 /pi, (rather than Mach number) in figure 1. A close empirical fit to the 
curves is given by 


where 


G 


= 3(y+l) 
10 
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Symmetrical Shock Solutions 



The shock solution is based, on flow con- 
tinuity considerations applied to the region 
between the body and its associated shock 
depicted in sketch (a). The length rep- 

resents the shock standoff distance from the 
body in the symmetry plane and, for a given 
free-stream flow, is related to the body and 
shock radii, and R g , as will be shown. 

The length A* represents the shock 
standoff distance from the sonic point on the 
body located at y*. The ratio A*/y* 
depends on the sonic point inclination angle, 
<f>*, for a given free-stream flow. This depen- 
dence establishes a unique relationship between <}>* and the shock inclination 
angle, 0*, as will be shown. 


As a consequence of the foregoing, if <j>* and y* can be specified for a 
body in a given flow, the parameters A s ^., R s , A*, and 0* are determined and 
constitute the shock solution. For this purpose, explicit relationships 
between the shock and body parameters are developed in the following sections. 


Shock shape — Reference 1 defined shock traces as conic sections of suit- 
able bluntness B s (see Notation) and obtained the shock solutions for a range 
of conic-section blunt bodies. In the present investigation, it was assumed 
that conic sections can closely approximate the shock trace for a larger vari- 
ety of blunt bodies than were considered in reference 1. The shock bluntness , 
B s , is determined from the shock parameters R s , y*, and 0* by 


B 

s 



( 5 ) 


Shock and body parameters at symmetry plane — A mass-flow continuity analy- 
sis applied to a "volume" element adjacent to the plane of symmetry in two- 
dimensional flow results in the following relationship for the shock and body 
radii, R s and R^, the shock standoff distance, A st ’ and the normal shock- 
density ratio: 



A detailed derivation of equation (6) is presented in appendix A. The values 
of A s ^./R s as a function of Ag^/R^ for various ratios P 2 /P 1 are presented 
in figure 2. 


3 



Shock and body parameters at sonic "plane " — The average mass-flow compo- 
nent normal to the control surface A* (sketch (a)) is denoted as (pv)*. Mass- 
flow continuity fixes the ratio A*/y* as follows: 


A* P 1 V 1 

y (pV)* 

n 


(T) 


The mean normal mass-flow component, (pV)*, was determined as follows: 
The value, pV, has the maximum value of 

P*V* = p 1 V 1 (p 2 /pi) (2/[ (y+l) - (y-l)(p 1 /p 2 ) ]W Y 1 


at the sonic point and its derivative taken in any direction is zero. It was 
assumed that the second derivative at the sonic point is also zero. The normal 
(vertical) components of pV at the shock and the hody are, respectively, 
p 1 V 1 [ ( Pq/pi )-l]sin 0* cos 0* and p*V* cos <j>*. These specifications on pV are 
accommodated by a four-term polynomial in x/A* which is integrated to obtain 
the normalized mean mass-flow value: 



Equation (8) indicates that, while primarily sensitive to the sonic angle 
((>*, the value for A*/y* is determinate only if a functional relationship 
exists between the sonic angle <}>* and the shock angle . 0*. This functional 
relationship is established by shock solutions for a family of widely varying 
body shapes, including a flat plate, a circularly rounded slab, and a circular 
cylinder. The shock solutions for a flat plate and circularly blunted slabs 
are discussed in the following sections. These solutions, along with the cir- 
cular cylinder solutions of reference 1, were then utilized to establish A*/y* 
and 0* as functions of <J>*. 


Fiat plate — In the case of a flat plate supporting a shock of bluntness 
B s , the sonic point is assumed to reside on the windward side of the plate 
corner, whence <J>* = 0. Simple relationships between B s , A s ^./R s , A*/y*, and 
0* with <j>* = 0 were proposed in reference 3 for an axisymmetric disc on 
P 

purely geometric grounds. These relationships also extend to two-dimensional 
flat plates and can be expressed in the present notation as follows: 


B 

s 


= V A st 


(9a) 


( A*/y*) = (A /R ) /tan 0* 

p st s p 

Equations (8) and (9b) are combined to give 


(9b) 


h 


4 tan 0* 

E 


( 10 ) 



(— - l) sin e* cos e* + 3 £ ~- 
\p x ) P P PlVi 


Since 4 t /R is known from equation (6) and since A^/R^ = 0, solutions for 
0p (eq. (10)7 and for (A*/y*)p (eq. (9b)) may be calculated for a flat plate. 

The resuits of such calculations are presented in figure 3 for (A*/y*)p and in 

figure 4 for 0* as functions of the normal shock-density ratio, P 2 /P 1 ; the 

specific heat ratio, y, is used as a parameter. 


The equation for a shock with bluntness B g can be expressed as 



tan 0* 

SB s tan 2 0* + 1 


( 11 ) 


In the case of a flat plate, y* = b and 0* = 0*. Equations (9a) and (ll) are 

combined to relate the centerline shock standoff distance, A ^ , to the plate 
semiwidth, b : ■ 



/B tan 2 0* + 1 

_! E 

B tan 0* 
s p 


( 12 ) 


The shock solution for a flat plate is now completely specified since B g , 

A gt /b, (A */y*) , and 0* are in hand. 

P P 

Circularly blunted slab — It is assumed that a circularly blunted slab can 
be specified so that at some free-stream Mach number, a circular shock hawing 
nearly the same radius as that of the body will result. In this circumstance, 
A* = A and 0* = <j>* , and with y* = R g sin 0* = sin 0*, equation (8) gives 



R 

s 



4 tan 0* 


1 sin 0* + 


p*V* 

Pl v l 


(13) 


The expression A s ^/R g = A ^/R^ is determined using equation (6) and has a 
unique value for a given free-stream flow. Thus, 0* = <j>* is readily deter- 
mined using equation (13) and the corresponding value of A*/y* is given by 


A * 

y* 


R 

s 


St 


sin 0* 


(14) 


Sonic angle , <j>* — In the foregoing, the shock solutions for a flat plate 
and a circularly blunted slab were determined once 4>* was fixed, that is, 
with <j>* = 0 and with <j>* = 0*. Shock solutions for other body shapes require 
general rules for selecting the appropriate sonic-point angle. The rules for 
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determining the sonic-point location and inclination angle follow those of 
reference 3. The arguments and pertinent equations are presented in appendix 
B. The essential results are given in figures 5 and 6. Figure 5 gives the 
reference sonic-point angle, 4>*, on a circular cylinder as a function of the 
reciprocal shock-density ratio, Pi/P 2 » using the specific heat ratio, y, as a 
parameter. Figure 6 gives the sonic-point angle ratio in the form 
sin <f>*/sin <j>* , either as a function of the body ellipticity ratio a/b or as a 
function of the ratio r/b (for a round-corner slab section). 

Relationship of A */y* and 0* to <}>* — A preliminary plot of the calculated 

A*/y* values with sin <p* for a flat plate, a circularly blunted slab, and for 

a circular cylinder determined from reference 1 suggested a uniform exponential 

increase of A*/y* with sin <f>*. Accordingly, the three available A*/y* 

values were utilized to correlate A*/y* with sin <j>* over the range, 

0 < d>* < d>*: 

T c 



The values (A*/y*) and (A*/y*) correspond to the circular cylinder and flat 

P 

plate, respectively. The value of sin 4>* is normalized by the circular 
cylinder value, sin (J>£. The value n was determined by substituting the 
values of A*/y* and <f>* from the circularly blunted slab solution into equa- 
tion (15). The results of evaluating the expression A*/y* - (A*/y*)^ as a 

function of sin <|>*/sin are presented in figure 7 with the free-stream con- 
ditions expressed in terms of the normal shock-density ratio, P 2 /P 1 . 

A preliminary plot of the three available 0* values with sin <f>* sug- 
gested a correlation curve of some complexity. Therefore, additional values of 
4)* were required to define the relationship between 0* and <f>*. Such values 
were calculated by utilizing the A*/y* - <p* relationship (fig. j) with equa- 
tion (8). The results are presented in figure 8, in which the term 0* - 0* 

is expressed as a function of the normalized sonic-point angle variable, 
sin 4>*/sin <j>£, with normal shock-density ratio, p 2 /pi, as a parameter. 

It is interesting that the correlations of A*/y* and 0* with <f>* in the 
forms presented in figures 7 and 8 are independent of the specific heat ratio, 
y • The application of these correlations to shock solutions in general will be 
demonstrated in a following section. 


Asymmetrical Flow Shock Solutions 

The pattern of asymmetrical shock is assumed to consist of an upper and a 
lower conic section of different bluntness, B s ; both sections, however, are 
assumed to have a common vertical tangent and radius of curvature, R s , at the 
shock apex as shown in sketch (b). The appropriate sections constituting the 
shock are found by iterative calculations so that, ultimately, both upper and 
lower shock elements satisfy A*/y* values based on mass-flow continuity at 


6 





their respective body sonic-point ordinates. In this regard, it is assumed 
that equation (6) is applicable although the stagnation streamline is curved. 
The curvature of the stagnation streamline, however, does introduce certain 
stream displacements denoted by y s and y £ in sketch (b). An approximate 
analysis of these displacements is presented in appendix C. The results of the 
analysis are presented in figure 9. 

The asymmetrical solution is applied to a class of two-dimensional bodies 
■typified by the section shown in sketch (c). The radii r u and r^ are the 
leading and trailing edge radii of the section whose chord length is c. These 
radii are tangent to a body arc of radius R^. The tangency points fix the 
indicated angles, s u and e^. The angle of attack, a, is measured between the 
chord and the free-stream direction. 

The technique of the asymmetrical shock solution is best demonstrated by 
the following calculative example with necessary equations introduced as 
required. 

Calculative example — The solution is applied to a wing-like body such as 
shown in sketch (c). The coordinate system most convenient to the solution is 
shown in figure 10 along with specific geometric details of the body in 
question. 

The solution is based on an initially assumed value for the span b, which 
locates the shock apex (fig. 10). An approximate value for b is based on the 
following: If the angle of attack a is 90°, the flow will be nearly symmet- 
rical and the shock apex will be located approximately at midchord, or, 
b = c/2. If the angle of attack is sufficiently small that a sonic point will 
reside at the point where r u and R are tangent , the circular cylinder shock 
solution applies with <f>* = = 90° - a - e u . The shock apex will then be 

located at b = r u . A linear variation of b with an angle of attack a is 
assumed with the above limits, giving the following result: 


b _ 

1 

1 

r 

u 

VO 

o 

0 

R 1 

c 

2 

_2 

c 

[e + 

L u c- 1 
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A trial shock solution can be calculated using the values c = 1.0, r u = 0.10, 
r £ = 0.05, \ = 1.92, e u = 10°, e £ = l6.5°,.and a = 65° (fig. 10). 

A reference ordinate is defined as 


Y = (c - r - r„)sin a + r = (l.O - 0.10 - 0.05)0.906 + 0.10 = 0.870 
m u £ u 

(17) 


The free-stream flow is ideal air at Mj = 00 (y = 1.4, P 2 /P 1 = 6.0 ). The 

basic shock parameters associated with the flow are as follows: 

(A*/y*) = 0.608 (fig. 3).; 6* = 19.45° (fig. 4); and <J>* = 44.5° (fig. 5). 

p P c 

The approximate shock tangency point (apex), b, given by equation (l6) can 
then be evaluated: 


or 


1 0.10 

90° - 65° 

2 “ 1.0 

10 ° + 44 . 5 ° 


0.317 


b = 0.317 (c = 1.0) 

The upper sonic angle, <j>*, is determined as follows: The ratio r u /h = 

0.10/0.317 = 0.316 takes the role of r/b in figure 6 where the ordinate value 
sin ij)*/sin cf>£ = 0.77 is found. The value for <j>* = sin -1 (0.77 sin 44.5°) = 
32.7° is calculated. This angle is compared with the inclination of the upper 
corner tangency point, ^ = a + e - 90° = 65 + 10° - 90° = -15°. The sonic- 
point angle criteria defined in appendix B require the larger of these angles 
to represent 4>*; therefore, <t>* = 32.7°. 

The upper sonic point coordinates are given by the following: 

Y* = Y - r (1 - sin <f>*) = 0.870 - 0.10(l - 0.540) = 0.824 (l8) 

u m u u 

X* = (c - r - rjcos a + r cos <\>* = 0.85(0.422) + 0.10(0.842) = 0.443 
u u £ u u 

(19) 


With the angle <(>$ specified, the ratio (A*/y*) and the shock angle 0* 

are determined as follows: The values (A*/y*) - (A*/y*) and 0* - 0*, asso- 

* P 


dated with <j>* = 32.7° (at P 2 /P 1 = 6.0), are found from figures 7 and 8 which 

give (A*/y*) - (A*/y*) = 0.109 and 0* - 0* = -2.60°, whence 

P P 


( A* /y * ) = 0.109 + 0.608 = 0.717 

and u 

6* = - 2 . 60 ° + 19.45° = 16 . 85 ° 


The lower sonic angle, is determined as follows: The lower span, 

opposing the upper span b, has the value, Y m +r £ -b = 0.870+0.05-0.317 = 0 . 603 . 
The ratio, r^/(Y m +r^-b) = 0.05/0.603 = 0 . 083 , takes the role of r/b in 
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figure 6 where the ordinate, sin <f>*/sin <J>* = 0.52, is found. The value for 
<f>f = sin" 1 (0.52 sin 44.5° ) = 21.4° is calculated. This angle is compared with 
the inclination of the lower corner tangency point , 

<f>^ = 90° + e £ - a = 90° + 16.5° - 65° = 4l.5° 

The previously mentioned sonic-angle criteria 1 requires the larger of these 
angles to represent <j>*; therefore, <f>* = 4l.5°. 

The lower sonic-point coordinates are as follows: 


Y* = -r sin <)>* = -0.05(0.663) = -0.033 (20) 

X* = r £ cos <f* = 0.05(0.749) = 0.037 (2l) 


The ratio (A*/y*)^ and the shock angle 0| will now he determined. The 

values (A*/y*) - (A*/y*) and 0* - 0*, associated with if)* = 41.5° 

P P 

( P 2 /p 1 = 6.0), can he determined from figures 7 and 8, which give 

( A*/y* ) - ( A*/y* ) = O.I 65 and 0* _ 0* = -1.30°; therefore, 

P P 


and 


( A*/y*) = 0.165 + 0.608 = 0.773 

0* = -1.30° + 19.45° = 18.15° 


Sufficient information is now available to estimate the coordinates X s 
sind Y s , of the shock "apex" 


X 


s 


X* + x* + A* = X* + y*(x*/y*) 
u u u u J u J u 


+ y*(A*/y*) u 


( 22 ) 


The value 


Y = Y - b 
s m 

y* = h - (Y - Y*) = 0.317 - (0.870 - 0.824) = 0.271 
u mu 


(23) 

(24) 


is the vertical displacement between the upper sonic point on the body and the 
shock apex. The value x^ is the horizontal displacement on the shock from 
the apex to the point where the shock inclination is 0$. The displacement 
ratio, (x*/y*) , is approximated by the "mean" shock slope, tan(0*/2). Equa- 
tions (22) and (23) can now be evaluated: 

X = 0.443 + 0.271 tan — + (0.27l) (0.717 ) = 0.677 

S t 

Y = 0.870 - 0.317 = 0.553 
s 

Occasionally , the selection rule may give an angle <j>$ > <f>*, in which 
case is assigned the value of <f>£. In this event, the location on the 

body where its surface inclination is <j>* fixes the lower sonic-point 
coordinates . 
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The approximate shock apex is located by the above coordinates on figure 10 and 
is denoted by the circular symbol. 

Attention is now given to defining the stagnation streamline. For this 
purpose, the length Aq in the free-stream direction between the shock apex 
and the body, the associated body surface inclination 4> q , and the radius 
at the body terminus of Aq are required. These values may be calculated with 
the geometry of simple analytic body sections but, in general, are most expedi- 
tiously scaled from an accurate drawing of the section. In this case, the mea- 
surements from figure 10 are Aq = 0.277 and <J> o = 21°; has the given value 
of 1.92 units. A simplified stagnation streamline consisting of two straight 
lines of equal length, one drawn horizontally and one drawn at the inclination 
4> o , respectively, are located as indicated by the dashed line elements in fig- 
ure 10. While this approximation lacks the qualitative feature of the actual 
uniformly curved streamline, its length, 

A s A 0 ( 1 + cos <f> o ) /2 = 0.277(1 + 0.934)/2 = 0.268 (25) 

S L 

is a good approximation to that of the actual stagnation streamline. The angle 
4> s -t at the point of the stagnation streamline impingement on the body is (in 
degrees ) as follows : 

A q sin 4> 0 

4> st * *0 - 57.3— — = 21° - 57.3(0.277/2) (0.357)/l.92 = 19.5° 

(26) 

The stagnation streamline value, Ag^/R^ = 0.268/1.92 = 0.140, and the corre- 
sponding value, Ag^/Rg = 0.200, is found from figure 2 with P 2 /P 1 = 6.0. The 
value Ag^/Rg = 0.200 with P 2 /P 1 = 6.0 gives the ordinate y s /A s t s: '- n < t > st = 
0.18 from figure 9. The vertical displacement y s of the dividing (stagna- 
tion) streamline from the shock apex is evaluated as follows: 

y = (y /A sin <j> )A sin <J> = 0.l8(0.268)sin 19.5° = 0.016 

s s SX SX S X SX f \ 


The effect of the displacement y s is to reduce the mass flow intercepted by 
the upper shock and to increase the mass flow intercepted by the lower shock. 
These mass-flow continuity considerations then affect the values of A* as 
follows : 

A* = fyj - y s )(A/y)* 
a* = (y* + y s )(A/y)» 

With the value of R s = A s -f- /( A s ^./R s ) = 0.268/0.200 = 1.340, an improved 
determination of the ratio x*/y* is made with the conic-section relationship: 


(28) 


x*\ tan 9* 

y*J 1 + (R s /y*)tan 0* 


(29) 


A check of the shock apex or tangency point location X g as given by both 
the upper and lower shock geometry is then made. For this purpose, the values 
of y*, (x*/y*) u , and (x*/y*)^ need to be established first: 
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( 30 ) 


y* = Y - b - Y* = 0.870 - 0.317 - (-0.033) = 0.586 

x, m x. 


With equation (29): 


tan 16.85° 


u 1 + q* 2 j ± tan 16 . 85 ° 


= 0.121 


£ 1 + 


tan 18 , 15 ° 
1.3U0 , 


= 0.187 


tan 18.15° 


The shock apex location X s is given for the upper and lower shock hy 

X = X* + (y* - y + y*fe) 

su u VJ u \y*J u u [y*/ u 

= 0.443 + (0.271 - 0.016)0.717 + 0 . 271 ( 0 . 121 ) = 0.659 


x st * *{ + + MfJ, + 

= 0.037 + ( 0.586 + 0.016)0.773 + 0 . 586 ( 0 . 187 ) = 0.612 


A discrepancy or mismatch of x S £~^su = 0*^12 - 0.659 = -0.047 is indicated. 
Since the above values for X g are primarily affected by and y$ , an incre- 
mental correction term, V, is applied to these values so that X g £ = X gu . 
Equations (3l) are modified accordingly and equated: 


x * + ( y # + 7) 1 - 7 — + (y* + V) -jr 

u u \ v* + V/Vy*/ u \ V* 


Equation (32) is solved for V (neglecting nonlinear terms in V) with the 
result 


X - X 
s it su 


-0.047 

0.121 + 0.717 + 0.187 + 0.773 


= - 0.026 


The above correction, V = -0.026, when applied to b = 0.317, automatically 
modifies the values, y* and y*, to the required values. 

An iterative calculation is then made with the modified value, 
b = 0.317 - 0.026 = 0.291, replacing the initial value, b = 0.317. The calcu- 
lations are repeated with the exception that the initial value 
(x*/y*) u = tan( 0 */ 2 ) , used in equation ( 22 ), is replaced by the value (x*/y*) u 
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from equation (29). This substitution is clarified by examining the computing 
form presented in table 1. This form contains the foregoing calculations and 
presents the numerical results of iterations toward a convergent solution. The 
final results include the ordinate of the stagnation point location Y s -^ on 
the body and the shock bluntness parameters B su and B gJi given by equation (5). 


Pressure Distribution 


The pressure distribution is found by first determining the velocity dis- 
tribution V/V* = V(s), where s = s/s*, and applying the perfect gas 
relationship 

p =t l M v 7 - 

to give P/P £ as a function of s. This application is expedited through 
figure 11 which gives P/P + as a function of V/V* for various values of the 
specific heat ratio, y. 


Y-l 


(34) 


The method, in principle, is as follows: The velocity distribution func- 
tion is first determined for what is considered an analytic (smooth) body 
approximating the actual body shape. A correction term is then added if neces- 
sary to account for the effect of surface curvature discontinuities. The 
desired function V(s) is subject to^the constraints V(0) = 0 and V(l) = 1; 
the associated derivatives V(0) and v(l) are then evaluated at the stagnation 
and sonic points. 


The value of the stagnation point derivative 
terms in equation (3) with V 1 /V* = {p 2 /p\) 1 ^ 2 • 


V(0) 


is found by rearranging 


V(0) 


d( V/V*) 
d(s/s*) 



(35) 


The values Ag^/R^ and s*/A s ^. are provided by the previously determined shock 
solution for the body in question. 


The sonic point velocity derivative is related to the sonic point pressure 
derivative by differentiating equation (34) with respect to s/s* and evaluat- 
ing the result at V/V* = 1: 


d(P/P st ) _ yP*d( V/V* ) 

d(s/s*) P d(s/s*) 
st 


The sonic pressure derivative, to a good approximation, is -s*/r (ref. 3). 
Substituting this value for the left side of equation (36) gives the value of 
the sonic point velocity derivative as follows: 


V(l) 


d(V/V*) _ P st S * 
d(s/s*) yP*r 


(37) 
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The next subsection presents the derivation of velocity functions subject 
to the above constraints and applicable to smooth bodies (no surface curvature 
discontinuities). The following subsection presents an approximate method that 
accounts for the effect of surface curvature discontinuities. 


Veloaity distribution for a smooth body — In reference 5 it was pointed out 
that in the vicinity of a sharp sonic corner, the velocity in the present 
notation is 

V(s ) ~ 1 - /I - 5 (38) 


Equation (38) satisfies the requirements that V(0) = 0 and V(l) = 1, and gives 
V(l) = °° in agreement with equation (37) for a sharp sonic corner (r = 0). 
Equation (38) is generalized in a simple manner to accommodate V(l) ^ °°, but 
to retain the property, V(0) =0 and V(l) = 1: 


V(s) 


1 - /I - k]S 
1 - A - ki 


(39) 


The appropriate constant, , is determined by differentiating equation (39) 

and equating the result at s = 1 to V(l) (eq. (37)): 


ki . W(i) [v(i) - 11 , (U0) 

[1 - 2V( l) ] 2 

Equation (39) with ki satisfies the condition V(0), V(l), and V(l), but only 
in fortuitous cases gives the correct value for V(0). An independent velocity 
function is required to satisfy_the conditions V(0), V(l), and v(0). A mono- 
tonically increasing function, V(s), similar to that of equation (39) is chosen 

A - k 0 (l - 5) - /I - k 0 ~ 

V{s) = ~ZZ=I (*+!) 

where 1 - A - kQ 


4V(0)[V(0) - 1] 
[1 - 2V ( 0 ) ] 2 


Equation (4l) satisfies the "end" conditions V(0), V(0), and V(l). All "end" 
conditions are satisfied by a linear interpolation between the values given by 
equations (39) and (4l): 


V(s) = 


A - k 0 (l - s) - /l - k 0 


1 - 'A - k 0 


(l - s) + 


1 - A - k ! 1 


1 - /l - k x 


(42a) 


In cases where V(0) < 0.5, equation (42a) will not satisfy both V(0) and 
V(l) =1. A truncated series term is substituted for such cases: 
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V(s) = [V(0)s + (1 - V(0))s 2 ](l - s) + 


1 - /l - k 


is 


1 _ A - k! 


(42b) 


Equations (42a) and (42b) give a uniform velocity distribution and are presum- 
ably applicable to a smooth body of uniformly varying surface curvature. 


Effect of curvature discontinuities — In cases where the location of the 
stagnation point is Newtonian (most forward point on the body) or when the 
absolute value of the body surface inclination angle increases with distance 
from the stagnation point, the effect of surface curvature discontinuities may 
be estimated in the following manner. Newtonian theory suggests the relation- 
ship v = sin <J> (the lower case v is used to distinguish velocities from 
those given by equation (42)). The further approximation v s <j> (<f> in 
radians), is adopted to simplify the analysis. The body surface is assumed to 
consist of three tangent radii, Eg, r^ , and r. The values of v =.<j> associ- 
ated with s = 0, the tangent points , and at s = 1 are, respectively, (f> s t » 

A', » and 4>* as shown in sketch (d). 

Straight-line elements are drawn between adjacent 
<j> values resulting in a line function indicated 
by the solid line. The slopes of the individual 
elements are determined by the appropriate sur- 
face^radius , that is , where the element radius is 
Rg, V = s*/R b , etc. A smooth curve approximation 
to the line function is determined as follows : 

The midsection of the line function is extended 
to s = 0 and s = 1. The result is a linear base 
function, vg, to which smoothly varying incremen- 
tal functions Av' and Av” of the form Av' = 
C'(l-s) n and Av" = C"(s) n are added to match 
the line function ordinates and slopes at s = 0 
and s = 1: 



v = v B - Av r + Av" (43) 

The smooth Newtonian velocity distribution, given by equation (43), is analo- 
gous to the velocity distribution given by equation (42). The difference in 
velocity distribution between the line function and equation (43) represents an 
additive correction term accounting for surface curvature discontinuities. It 
is assumed that this Newtonian correction, suitably scaled, can be applied to 
the distribution given by equation (42) to account for the effect of surface 
curvature discontinuities . 


Detailed derivations giving the explicit form of equation (43) are omitted 
herein. Instead, the results of such derivations are embodied in the computing 
form presented in table 2. Additional instructions are given in appendix D. 

Numerical example — In graphical form, figure 12 shows sample numerical 
results obtained with equation (42) for the velocity distribution on a smooth 
body equivalent to a blunted cone. 2 The appropriate velocity correction term 

2While the shock solution presented herein applies only to two-dimensional 
flow, the method of determining the velocity distribution is general. The 
velocity gradients for this example (fig. 12) are determined from reference 3. 
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is provided by a calculated example in the computing form of table 2. Comments 
on the selection of suitable scaling factors and other details pertaining to 
the computing form are discussed in appendix D. 

The pressure distribution corresponding to the corrected velocity distri- 
bution was determined from figure 11 and is included in figure 12. 


COMPARISON OF EXPERIMENTAL AND PREDICTED RESULTS 


The validity of the present method was assessed on the basis of compari- 
sons with the theoretical shock envelopes of reference 2 and comparisons of 
pressure distributions with the theoretical results of references 1 and 5 and 
the experimental results of references 6 and 7. 


• Shock Shape 

; Comparisons of shock shapes predicted by the present method and reference 
2 for a flat plate are shown in figure 13. Good agreement is shown at the 
angles of attack, 90° and 76° (M} = 20). The agreement deteriorates somewhat 
at an angle of attack of 60°. A point to be noted at the angle of attack of 
60° with M^ = 5 is the location of the stagnation streamline. The streamline 
predicted by reference 2 is essentially straight and parallel to the free- 
stream direction, implying that the stream entered the shock at its apex. This 
is in contradiction to the assumption made herein that the stagnation stream- 
line should exhibit a subsonic type behavior and curve toward a nearly normal 
impingement to the body surface as indicated. Unfortunately, no experimental 
data are available to resolve the difference in these predicted streamlines. 


Surface Velocity and Pressure Distributions 

A comparison of the velocity distributions for a flat plate at Mj = 20 
predicted by the present method and those of reference 2 for a flat plate are 
shown in figure l4. Good agreement is shown between the predicted (positive) 
velocities leeward of the stagnation point at the angles of attack, 90° and 76°. 
There is less agreement between the predicted values when the angle of attack 
is 60°. The lack of agreement between the predicted values of the (negative) 
velocities windward of the stagnation point is more pronounced at the angles of 
attack, 7 6° and 60°. This lack of good agreement arises partly from the dis- 
crepancies between the predicted streamlines, since the windward velocity is 
very sensitive to the stagnation streamline as the stagnation point approaches 
the plate (windward) leading edge. 

Predicted and experimental pressure distributions over the forward face of 
various blunt bodies are compared in figure 15. The comparisons are made for 
axisymmetric bodies since two-dimensional data, particularly for asymmetrical 
shapes at various angles of attack, are relatively scarce. However, since the 
method of estimating the pressure distribution is independent of the 
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dimensionality of the flow, the good agreement shown in figure 15 between the 
results of the present method, the results presented in reference 5, and the 
experimental data demonstrates the validity of this method. 

CONCLUSIONS 


A method has been developed for predicting shock envelopes and pressure 
distributions for a variety of two-dimensional blunt bodies at various angles 
of attack. The method is restricted to those cases in which the bow shock is 
detached from the body and the flow over the forward face of the body is 
subsonic. 

This method utilizes a quickly convergent calculative procedure based on 
correlation functions which relate the shock standoff distances at the stagna- 
tion and sonic points to the body geometry. The correlation functions depend 
primarily on the normal-shock density ratio modified to a small degree by the 
specific heat ratio of the gas. Since the effect of the specific heat ratio is 
small, the method should give adequate solutions for the equilibrium flows of 
real gases. 

Predicted shock envelopes and pressure distributions compared favorably 
with the results of other theoretical methods and with experimental values in 
the Mach number range 4.0 to 20. 


Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, Calif. 9^035, Sept. 25, 1972 
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APPENDIX A 


BODY-SHOCK GEOMETRY IN SYMMETRY PLANE 


The relationship between the shock and body radii, R s and R-^ , and the 
shock standoff distance, A g ^. , in the plane of symmetry is determined by utiliz- 
ing mass-flow continuity, oblique shock relationships, and a suitable 
coordinate system. 



m 


Mass' flow near the shook — The 
shock and flow geometry is depicted in 
sketch (e). The analysis is confined 
to the domain below a control surface 
m-m which extends normal to the body 
surface (not shown). This surface is 
inclined to the symmetry plane at a 
small angle, <f> . The purpose of the 
following calculations is to find the 
mass-flow component normal to the plane 
m-m and its derivative along m-m at 
the shock. 


A streamline entering the curved 
shock of radius, R s , at some point, y^, 
above the shock apex will be deflected an angle, 6^, and will pursue a curved 
path of radius, <*5. , behind the shock. The results of reference 8 give the fol- 
lowing radius ratio for small 0: 


R 

s 
c R \ 


P 2 
Pi 


P 2 


■ , (5 ~ y) 

Lp ! + y)_ 


0 


(Al) 


where 0 = y/R s and p 2 /Pi is the normal-shock density ratio. 


By oblique shock theory, the small deflection angle 6^, can be expressed 
as follows: 


= ( W 


/P2 

\p7 


- 1 


(A2) 


In the vicinity of the shock, the equation of the streamline entering the 
shock (x = 0) at y = y^ is given by Maclaurin's expansion 


y = 


n + if-ir 

l Pl / R C 


P 2 

CM 

P 

2P1 

Lpi 


+ (5 - y) 

(l + y) 




(A3) 


The equation of the control surface, m-m, is 


y = y s - 4>x 


(AU) 
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The inclination of the streamline where it intersects the control surface 
is dy/dx + <f> . Flow continuity along the streamline, s = x, gives the mass- 
flow relationship, pV/piVi = ^/y. The component of mass flow normal to m-m 
is 



(A5) 


Equation (A5) is expressed in terms of y g and x as 


■O 

<!_ 

S3 

1 

H 

/y s - x 

(fO 

\Pl 

p 2 

1 

P 2 , (5 - yfl 

' 

X 

PlVj f(x)\ 

> f(x) 

R 

s 

Pi 

_Pl (1 + y)_ 

B S 2 ] 


+ <j> 


Equation (a 6) gives at the shock (x = 0) 

(pV) 


n 


= 6 + <i> 

PlVi s 

The derivative of equation (A 6) at x = 0 is 

d(pV) 


L(pV) /p \ 

= - r[—~ 7 (6 s + + ^r 


(A 6) 


(AT) 


(a8) 


A two-term expression for the normal mass-flow component in the vicinity 
of the shock is 


Pl»i 


(s + <f>) - 

s 


2/ P 2 

LM P1 



+ <J>) 



(A9) 



Mass flow near the body — A streamline in the 
vicinity of the stagnation point on the body is 
shown in sketch (f). At a small distance, dx, from 
the body, the stream is essentially parallel to the 
surface. 


From equation (3), the mass-flow gradient at 
the stagnation point of the body is 


d(pv) n a * A st /y 

PI 7 ! 45 " 4 st 


(A10) 


The subscript n has been appended to indicate that the flow at (or near) the 
body surface is normal to the control surface m-m. 


At the point ds = R^ on the body surface, the mass flow is 
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(All) 


(pV) (1 + A ,/R, ) 

-^7 ^ °V 

Assuming that the circulation <j> V • ds = 0 holds for the perimeter of 
small concentric strip of width dx and length R^<J> it can be shown -that 

d(pV) (l + A./R) 

-—a -si-i- g* ti 

p^dx A 


Normal mass- flow variation along the control surface — A four-term poly- 
nomial expression for the variation of normal mass flow from the shock to the 
body along the control surface m-m is 

s - [•& - *) * " x ■ • ‘-fey * "•(*')■ 


The constants C2 and C3 are determined so that x/A gt = 1 (body surface) , 
equation (A13) reduces to equation (All) and the derivative of equation (A13) 
reduces to equation (A12) . 

Shock and body radius relationship — The normal mass-flow distribution 
given by equation (A13) was integrated over the interval 0 < x < to give 

the total mass flow passing through the control surface m-m and was equated 
to the flow entering the shock element, y s PiVj. Omitting arithmetic details, 
the following identity resulted: 


rU +^r -i/P? \ A 


( A sty ^ x\Gi 1 4 st 


The geometric relationship y s /R g = (Ag^+R^ )<J>/R S relates 6 g , ()>, and 6 S . The 
ratio A s ^./^ and the angles 6 g and <j> were thereby eliminated from equation 
(Al4) , giving the final result: 


- 1 + 


X ) 1 - 




3 VW 


jt If l ( Q J\ 2 

R s ~ 2 v v 3 W 


xtiri + 2 
Y + l\Pl/ 3 
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APPENDIX B 


BODY INCLINATION AT THE SONIC POINT 


The solution for an infinite shock-density ratio was used to determine the 
inclination of the sonic point on a conic-section body. The results were then 
modified to apply to finite shock-density ratios and to nonconical sections by 
the following procedures . 


Conic-Section Body 

The inclination angle <J>* at the sonic point was determined as a function 
of the body bluntness parameter B^ using the Busemann solution for an infi- 
nite shock-density ratio. Details for the axisymmetric flow solution are pre- 
sented in appendix A of reference 3. A similar solution for two-dimensional 
flow gives the following result: 

(l - sin 2 <(>*) - x - y { [ (B^ - l)-sin 2 <)>* + l] 3 / 2 - [(B^ - l)sin 2 <|)* + l]) 

(Bl) 

In the present application of equation (Bl) , the gas was assumed to be 
perfect. Since the theory requires an infinite shock-density ratio, the value 
P*/P g -fc = 0.607, appropriate for a gas with y = 1.0 and Mi = °°, was used. The 
results for sin <j>* are normalized by the value for a circular cylinder 
(sin <|>* at B^ = l) and are presented in figure 6 as a function of a/b (note 
b/a = B^/ 2 ). Although the absolute values for sin <j>* given by equation (Bl) 

may not be correct for a finite shock-density ratio, it is assumed that the 
ratio sin <f>*/sin <j>£ will not differ greatly from that given by exact theory. 
Thus, if the value of <f>* is known for a circular cylinder in a given free- 
stream flow, the value of 4>* for a conic-section body can be determined from 
figure 6. The values of <f>* for a circular cylinder by the exact solutions of 
reference 1 are presented in figure 5. 


Nonconic-Section Bodies 

If a blunt body does not have a conic section, the parameter B^ is not 
appropriate. However, the sonic-point location on other types of blunt bodies 
may be estimated on the basis of the foregoing analysis. To this end, it is 
assumed that a shock in the form of a conic section is appropriate for all 
classes of blunt bodies considered here. For example, a blunt body with a 
rounded corner of radius r is illustrated in sketch (g). When the flow has 
an infinite shock-density ratio, the shock coincides with the body in the 
vicinity of the sonic point. If the sonic point is on the rounded corner of 
the body, the Busemann solution applies if the body corner radius r is 
equated to the radius of a conic section at the same sonic-point ordinate, y* 
(ref. 3): 
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r_ 

b 


(B2) 



Sketch (g) 


1 

1 + (B^-l)sin 3 <{)* 

Equation (B2) gives the r/b value for an 
equivalent body that supports the same elliptical 
shock and the same sonic angle as a body with the 
bluntness parameter B^. The value of 
sin <|>*/sin <j>£ for rounded corner bodies as a 
function of r/b is plotted in figure 6. 

In some cases, a qualification of the sonic 
angle is necessary, for example, when the body is 
a large-angle wedge or a large-radius surface 
terminated by a rounded corner. In such cases, 
the inclined or circularly blunted portion of the 
body will be tangent to the corner radius r at 
which point a surface inclination angle is 

defined. Tfye sonic-point angle is then taken as 
(jjt", if <j)t i s greater than the value of <f>* 
determined from figure 6. 
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APPENDIX C 


ASYMMETRIC STAGNATION STREAM 


Streamline Impingement Angle on Body- 


Consideration is given to a small domain in the vicinity of the stagnation 
point on the body (sketch (h)). Because of vorticity effects induced by the 
shock, the stagnation streamline does not impinge normally to the body but will 

be inclined at some angle ft with respect to the 
normal. In this small domain the flow is incompres- 
sible and the velocity components will be linearly 
proportional to the local coordinates , 5 and n . A 
stream function which yields linear velocity compo- 
nents and satisfies the boundary condition \p = 0 
at 5 = ftn and n = 0 is 



^ = c(5n - ftn 2 ) 


(Cl) 


The velocity gradient 9u/95 was assumed constant 
along the stagnation streamline of length A s -j. , and 
to the order of accuracy required here, was 
assigned the value 


3 2 ij> _ 3_u 

953n 95 


— VxA , = C 
P2 st 


(C2) 


The vorticity 5 (ref. 9) in terms of the coordinate system used herein is 

$r-& • ‘ - T >(£n i - - <«> 


where 0 is the small angle indicated in sketch (e). 
s 

Equation (Cl) is differentiated partially for 5, and with the C from 
(C2), gives 


~ If = frf- + 777 = - 2C ,,, n = 2 Pl Vift/p 2 A 


9u 


9n 95 


35 




st 


(CU) 


Equations (C3) and (C4) are equated to yield 


2 ^\pi 7 E s 


(C5) 
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Stream Geometry 

The stagnation streamline is referenced to a 
coordinate system whose origin is at the shock 
"apex" as shown in sketch (i). The streamline ordi- 
nate is y s at the shock and y s ^. at the body. 

Examination of incompressible asymmetric flows 
in general shows that the stagnation streamline 
increases uniformly in curvature as it approaches 
the stagnation point. A simple curve having this 
property was chosen to represent the stagnation 
streamline between the shock and the body. 


^ = C^/tan 6 (C6) 

The constant C is selected to give the correct streamline radius at the 
shock as follows : 

C 'yf = ^ tan 6 = (C7) 

<*• s s 

The differential streamline length is 


^ dlS 

tan 6 


( C8) 


After integrating equation (c8) and evaluating it between the limits 6 S 
and <t> 3 ^. + ft, it is equated to the length of the stagnation streamline, A st , to 
yield 

A . = $ tan 6 Jin 
st s s 


smif , + u) 

st 

sin 6 


(C9) 


Similarly, the ordinate change in the streamline between the shock and the body 
can be found by integrating dy = sin { ds = cos 6 d6: 


y , - y = R tan 6 [sin(4> . + ft) - sin 6 ] 

SC S S S S u S 


(CIO) 


The term R s can be eliminated between equations ( C9 ) and (CIO) and, using the 
relationship = [ ( P 2 /p l )-l](y s /R s ) , the following streamline ordinate 
relationship results : 


y st - y s 

A + 
st 


sin(4 l . + ft) - sin 6 
st s 


£n 


sin(<|> , + ft) 

st 

sin 6 


(Cll) 


Equations (C5) and (Cll) were used to calculate values of y s and y s ^. as func- 
tions of A s .j. /R s , (}i s ^ , and P 2 /P 1 . The results are presented in figure 9. 
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APPENDIX D 


AUXILIARY INSTRUCTIONS FOR TABLE 2 


The computing form (table 2) for determining the velocity correction term 
applies directly to a body consisting of three finite tangent radii, R^, r^, 
and r (items (l) to (3)) between the stagnation and sonic point locations 

(sketch (j)). The angles t ^ 

(items ( k ) to (7)) measure in radians the surface 
inclinations at the stagnation point, the tangent 
points, and the sonic point with respect to a plane 
normal to the free-stream direction. In cases 
where any of the surface element radii are of zero 
curvature, the distances , s^ , and s* to the 
tangent points (items (8) to (lO)) as specified in 
the computing form are indeterminate. In this 
event , the following auxiliary equations for com- 
puting the values of s are required. For clari- 
fication, the geometry for each of the following 
cases is shown in sketch (k). 



Case I, R^ = °°, r and r Finite 

Since the angle 4> g ^ = ^ , the value of (item (8)) is 

t" t" 

. , r^sin <j> - sin <t> , ) - r(sin (|>* - sin <J> ) 

"C «. St 

s = y* - ■ — ■ . .. ■ ■ - - ■ 1 

cos 6 . 
st 


(Dl) 


Case II, r = °°, R^ and r Finite 
t* t" t" 

Since the angle 4> = <j> , the value of s (item (9)) is 


t I ■£ t 

7* - R t) (sin $ - sin (f> gt ) - r(sin <j>* - sin 4> ) 

s = s + £7 

cos <f) 


(D2) 


Case III, r = °°, R^ and r Finite 
t" 

Since the angle <p = <p*, the value of s* (item (10)) is 

y* - R b (sin - sin - r^(sin 4>* - sin <t>^ ) 


5 * = 


COS <p* 


(D3) 


2k 




Case IV, R^ Finite, r 1 = r = » 

t 1 t" t" 

Since the angles <j> = <j> = 4>* and the value of s = s* (item (9) 

equals item (lO)), 


t" * f y* - \(sin - sin 4 > st ) 
s = s* = s + 


cos <f>* 


(DU) 


Case V, r"*" Finite, = r = 00 

t » t” t ' 

Since the angle <|> , = <j> and the angle <p = <j>* , the value of s 

(item (8)) is 

y* sin <J>* - Ax* cos <j»^ - r^[l - cos(<J>* - (j)^)] 

s sin(4>* - <(> ) 


(D5a) 


where Axg is the horizontal distance between the stagnation and the sonic 
point locations. 


The value of s* (item (10)) is 


11 y* - s cos (j> st - r (sin 4 >* - sin <|> o+ ) 


s* = s“ + 


st 


cos 4 )* 


(D5b) 


Since the angles 
equals item ( 9 )) is 


Case VI, R^ = 


t 

r 


r = Finite 


p st 


t ' t" 

4> = <)> , the value of 


t* 

s 


S t (item (8) 


t' 

s 


t" 


y* - r(sin <j>* - sin 


cos (j> 


st 


(D6) 


Scale Factor 

The "smoothed" Newtonian velocity distribution, equation (43), is not gen- 
erally equal to that of the velocity distribution given by equation (42). It 
is assumed, however, that the incremental velocity distribution, accounting for 
the effect of surface curvature discontinuities by Newtonian theory (item (25) )» 
is qualitatively valid and by an appropriate scaling factor can be applied as 
an additive correction term to the "smooth" velocity distribution given by 
equation (42). ^A, logical scaling. factor , f, is the ratio of the velocity 
gradients, f = v/v, where V and v are the respective gradients of equations 
(42) and (43). Preferably, the gradients should be evaluated at the 
stagnation point, whence 

f = ?(0 )/v(0) (D7) 
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The term V(0) is given by equation (35) and v(0) = s*/R^ is item (ll) in the 
computing form. 

In the case when R^ = 00 , the Newtonian value, v(0) = 0, is not realistic. 
An alternative choice is then the velocity gradient ratio at the sonic point: 


_ V(l) ? st 
' v(l) ’ ^ 


(D8) 
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TABLE 1.- SHOCK SOLUTION COMPUTING FORM 


c (chord) = 1.00 


r u = 0.10 


r £ = 0.05 


5 e £ = 16.5° 

6 a = 65 ° 


(©-©-©) cos © = 0.359 13 

Y m = ©tan©+© = 0.870 l!» 

= “ 15 

y = 1.1* l6 

[ 1+ ©]/[©-!- = 6 17 

G, fig. 1 = 0.966 18 



19 

Initial b = ©(l) 

0.317 

Zjb. 29 1*^0. 29 2 

1*8 

a o 1 

0.277 

0.2U0 


20 

r u /b 

= ©/© 

.316 

.31*1* 

1*9 

Scaled from fig,J0 
< * >0 at coordinates ©) - ©) 

21° 

20° 


21 

© 

ifig.6, 2I24I 

0 * sin <J>* 

.77 

.79 

50 


1.92 

1.92 


22 

23 

= 

<f>* = 

sin -1 ( ©) sin ©)) 
larger of (22) & (16) 

32.7 

32.7 

33.6 

33.6 

51 

52 


.268 

.11*0 

.233 

.121 


2h 

Y* = 
1 u 

(Jj) -(^23 ( 1-sin (23|) 

.821* 

.825 

53 

(52) &fig. 2. dgt/Rg 

.200 

.201 


25 

xa = 

©<2)cos © 

.1*1*3 

.1*1*2 

5l* 

26. 7[ ©sin©]/© 

1.1*° 

1.1° 


26 

sin 1 

© /sin © 

.77 

.79 

55 

♦st = ©-© 

19.6° 

18.9° 


27 

© 


.109 

.115 

56 

©) & ord. fig. 9(a) 

.180 

.178 


28 

© 

& fig. 8, e*-0« 

-2.60° 

-2.55° 

57 

y s = ©© sin © 

.016 

.013 


29 

£ - © + @ 

.717 

.723 

58 

Rs = ©/© 

1.31*0 

1.159 


30 

®5 = 

@ + © 

16.85° 

16.90° 

59 

l+[ © tan©]/© 

2.1*93 

2.1*31 


31 

r z /b 

= ©/©+©-©) 

.083 

.080 

60 

(p-) u =[tan©]/@ to©0 

.121 

.125 


32 


„ sin 

&fig.6, s . n ^ 

.520 

.510 

61 

y* = ©-© 

.586. 

.612 


33 

<(>* = 

sin -1 ((© sin(lj^) 

21.1* 

20.9 

62 

l+[@)tan©]/© 

1.750 

1.621 


3!* 

** = 

larger of ©) & @) 

1*1.5 

1*1.5 

63 

(pr)* = tan©/© 

.187 

.202 


35 

n = 

-©sin © 

-.033 

-.033 

61* 

A 5 = (© - © (§> 

.183 

.169 


36 

x* = 

(J) cos (© 

.037 

.037 

65 

A* = (©+©)© 

.1*65 

.1*83 


37 

sin 1 

© /sin © 

.91*5 

.91*5 

66 

© + ©© + © 

.659 

.61*2 


38 

© 

«is.T. (£) p 

.165 

.165 

67 

© + ©© + © 

.612 

.61*1* 


39 

© 

&fig.8, e*-e* 

-1.30° 

-1.30° 

68 

= © - © ' 

-.01*7 

.002 


1*0 

if) 

t = © + © 

.773 

.773 

69 

© + ©+©+© 

1.798 

1.823 


1*1 

9 * = 

©+© 

18.15° 

18.15° 

70 

7 = ©/© 

-.026 

.001 


1*2 

y* = 

© + ©-© 

.271 

.21*6 

71 

b = © + ©to©@ 

.291 

.292 


1*3 

*5 = 

@© 

.191* 

.178 

72 

(©and ord. fig. 9(b) 



0.680 

1*1* 

app. 

(?)u = tan( ® /2) 

.11*8 

0.121 0.125 

73 

Y st =©+©@sin© 



.630 

1*5 

*5 = 

@© 

.01*0 

.030 

71* 

B su = (@! /©) 2 =cot2@ 



11.36 

1*6 

Xs = 

©+© + © 

.677 

.650 

75 

B st = (@/@) 2 -cot 2 © 



-5.72 

1*7 

Y s = 

©-© 

.553 

.579 







29 








TABLE 2.- COMPUTING FORM FOR VELOCITY INCREMENT, AV 


R b =i.O 

I’"t'=oo 


4 4> st =0 a 

5 ♦ t *=0.524 

5 <j> t "=0.524 

7 cf)*=0. 524 


3 s t, = 0-0)© =0.524 

9 s t "=(8)+(©-0)@ =4.564 
3 s*=®+(0-@) 0=4.564 


0 

0.04 

.08 

/© =0.1146 
.15 

W -25 
/(®=1.0 
.60 
. .80 
1.0 




2 


Line 

function 


0 

0.157 

.281 



0 

0.026 

.085 



1 - (20 



0 

.04 

.08 

.1146 

.15 

.25 

1.00 

.60 

.80 

1.00 


1.000 

.96 

.92 

.8854 

.85 

.75 

0 

.40 

.20 

0 


1.00 

.701 

.484 

.346 

.243 

.092 

0 


0.52 
.367 
.253 
.181 
.127 
.048 

0 



0.52 

.367 

.253 

.181 

.127 

.048 

0 



a Note: 0, in radians. 

^Scaling factor f = 0.64, equation (D7) or (D8) 
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Figure 1.- The G function. 
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Figure 2.- Values of A st /R s as a function of Ag^/R^ for various values of Pz/Pl- 



(A*/y 



Pz/P\ 


Figure 3.- Sonic-point standoff distance (A*/y*)p for a flat plate normal 

to the stream. 
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Figure 4.- Shock inclination angle 


for a flat plate normal to stream. 



P/P 2 


Figure 5.- Sonic-point angle <f>* on circular cylinder. 



a/b, r/b 


Figure 6.- Sonic-point angle <p* on elliptic and round-corner slab sections. 
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.4 .6 

sin <p*/ sin <p * 


Figure 7.- Values of (A*/y*) - (A*/y*) as a function of sin 4>*/sin <j>£. 
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(a) y s 
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(b) y st 
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Figure 9.- Streamline displacement, y g and y s ^. 
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Figure 10.- Coordinate system for shock solution with example solution 
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Figure 11.- The ratio 



Figure 12.- Sample numerical results for a blunted cone, Mj = b.86. 





Figure l4.~ Comparison of predicted velocity distribution on a two-dimensional 

plate of various angles of attack. Mi =20. 
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s/s* 


Figure 15.- Predicted and experimental pressure distributions on various bodies 

NASA-Langley, 1973 — 1 A-^564 ^5 
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